Prediction of femtosecond oscillations in the transient current of a quantum dot in the 

Kondo regime 
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We invoke the time-dependent non-crossing approximation in order to study the effects of the 
density of states of gold contacts on the instantaneous conductance of a single electron transistor 
which is abruply moved into the Kondo regime by means of a gate voltage. For an asymmetrically 
coupled system, we observe that the instantaneous conductance in the Kondo timescale exhibits 
beating with distinct frequencies, which are proportional to the separation between the Fermi level 
and the sharp features in the density of states of gold. Increasing the ambient temperature or bias 
quenches the amplitude of the oscillations. We attribute the oscillations to interference between 
the emerging Kondo resonance and van-Hove singularities in the density of state. In addition, we 
propose an experimental realization of this model. 

PACS numbers: 72.15.Qm, 85.35.-p 
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There is a growing need to provide a quantitative de- 
scription of the sudden switching behaviour of single elec- 
tron transistors, since the continuous shrinking of conven- 
tional MOSFET transistors* and developments in state- 
of-the-art nanotechnology experiments suggest they may 
one day constitute the building blocks of organic comput- 
ers. Real-time electron dynamics in these devices have 
also profound implications for quantum computing 2 - and 
the realization of an electron analog of a single-photon 
gurA 

The effect of sudden perturbation in the form of step- 
like switching of the gate or bias voltage has been stud- 
ied thoroughly 4-7 and it has been unambiguously shown 
that the resulting transient current exhibits different time 
scales^—. The initial fast non-Kondo timescale is char- 
acterized by the reshaping of the broad Breit-Wigner 
resonance located around the impurity level, while the 
sharp Kondo resonance pinned to the Fermi level reaches 
a metastable state at the end of the much longer Kondo 
timescale. Subsequent investigations predicted that the 
asymmetric coupling of the dot to the contacts may in- 
duce interference between the Kondo resonance and the 
discontinuities in the density of states of the leada^. Ex- 
tension of the diagrammatic Monte Carlo method 1 ^ to 
impurities out of equilibrium 14 verified the transient cur- 
rent's dependency on the bandwidth of the contacts^. 

Ramifications of the band structure of contacts in time- 
dependent transport have been elucidated previously for 
prototypical systems by directly solving the Green's func- 
tions in the time domai n 12 ' 16 ' 17 . Ab-initio calculations 
also suggest that the crystallographic orientation of elec- 
trodes may influence transport properties 1 -^. 

The purpose of this letter is to provide a more realis- 
tic picture relevant to an actual experiment by using the 
real density of states of gold obtained from first principles 
calculations as input in a study of the transient current 
through a quantum dot asymmetrically coupled to gold 
leads. The dot level is abruptly switched such that the 
Kondo resonance is present in the final state. We analyze 



the effects of the density of states of the conduction elec- 
trons in the contacts, the position of the dot energy level, 
the asymmetries in the couplings, and the temperature 
on the instantaneous current. 

The physics of this system is captured sufficiently by a 
single impurity Anderson Hamiltonian which describes a 
single doubly degenerate level of energy dot attached to 
continuum electron baths. There is a one-to-one mapping 
between this model and the Kondo model via Schriefer- 
Wolff transformation for the parameter range we will be 
interested in this paper. The auxiliary boson transforma- 
tion is performed for the Anderson Hamiltonian. In this 
method, the electron operator on the impurity is replaced 
by a massless boson operator and a pseudofermion opera- 
tor. The U oo limit is achieved by restricting the sum 
of the number of bosons and pseudofermions to unity. 
The Hamiltonian is then converted into 



H(t) = ^\e k n kaa + V a (e ka ,t)c\ a(y tf f a + 

ka<r 



(1) 



where f\ (f a ) and c\ aa (c kaa ) with a — L,R create 
(annihilate) an electron of spin a in the dot and in the 
left (L) and right (R) gold leads, respectively. Moreover, 
n a and n kaa are the corresponding number operators, V a 
are the tunneling amplitudes for the left and the right 
leads, and b' (b) creates (annihilates) a massless boson 
in the impurity. 

Assuming that the hopping matrix elements have no 
explicit time dependence, the coupling of the quantum 
dot to the contacts can be parametrized as r^(m(e) = 
fi(fl)^i(H)(e), where T L{R) = 2ir\V L{R ){e f )\ 2 is a con- 
stant and £_L(.R)(e) is the density of states function. In 
order to describe the density of states of the gold con- 
tacts accurately, we employ ab-initio calculations. 

The ab-initio calculation of the gold density of states 
is based on density functional theory. We make use of 
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FIG. 1: The total density of states of gold is shown by a 
black dashed curve as a function of the separation from the 
Fermi level. The red curve corresponds to the best fit to the 
ab-inito data by a linear combination of Gaussian functions, 
which are represneted by green curves. 



the well-established WIEN2K package^ ., in which the 
full-potential linearized augmented plane wave method is 
implemented with a dual-basis set. Specifically, the gold 
core states are defined by the electronic configuration Kr 
4d 10 4/ 14 5s 2 , whereas the valence states comprise 5p, 5d, 
6s, and 6p orbitals. The face-centered cubic unit cell of 
gold (space group Fm3m, lattice constant a = 4.080 A) 
with 4 atoms per cell is used in the calculation. For the 
exchange-correlation potential the generalized gradient 
approximation within the Perdew, Burke, and Ernzerhof 
(GGA-PBE) parametrization is adopted 21 . The radius of 
the gold muffin-tin sphere is set to 2.5 Bohr radii. More- 
over, the plane wave cut-off for the scalar relativistic basis 
function is given by R m tK max = 8 and l max = 10. In- 
tegrations in the reciprocal space for self-consistent field 
cycles apply the tetrahedron method and 752 k-points in 
the irreducible wedge of the Brillouin zone. The resulting 
density of states (DOS) is shown in Fig. Q] 

In order to represent this complicated curve in a func- 
tional form, we use a fitting procedure with a linear com- 
bination of Gaussian functions given by 



where 



Pi(e) 



exp(-2(- 



-?)■ 



(2) 



(3) 



The outcome of the fitting procedure, that involves six 
different Gaussians of varying linewidth and peak posi- 
tion, is also shown in Fig. [1] The number of Gaussians 
corresponds to the number of gold states in the vicin- 
ity of the Fermi level (five d orbitals and one s orbital). 
The resulting best fit parameters are summarized in Ta- 
ble |U Throughout the rest of this letter, we will switch 



to atomic units, where Ti 
calculations accordingly. 



e = 1, and perform our 



TABLE I: Fitting parameters for 6 different Gaussians used 
to replicate the density of states of gold (in eV). 
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We invoke the well tested non-crossing approximation 
(NCA) to obtain the pseudofermion and slave boson self- 
energies. NCA gives reliable results for dynamical quan- 
tities except for very low temperatures or finite magnetic 
field. We stay away from both these regimes. We solve 
the resulting real-time coupled intcgro-differential Dyson 
equations for the retarded and lesser Green's functions in 
a discrete two-dimensional grid. A technical description 
of our implementation has been published elsewher e 11 ' 22 . 

The net current in the electrical circuit can be derived 
from the Green's functions Gp S { H\t 7 t') and B< R \t,t'). 
We will denote the net current by J(t) = — Iil(t), 

where Ih{t) (/b(*)) represents the net current from the 
left (right) contact through the left (right) barrier to the 
dot. The general expression for the net current 2 ^ can be 
rewritten using pseudofermion and slave boson Green's 
functions^ 2 -, finally leading to 



I(t) = -2(T L - T R )Re [ I dhUt, h)h(t - h) ) + 



2T L Re (J <fti(&(t,ti) + &(t,ti))/L(t-fi) 
2f R Re ( f dh(Ut, ti) + £„(«, h))f R (t - h) 



(4) 

withf (t,ti) = G^fttO-B^ti.t) as well as&,(t,*i) = 
G* eu (t,t 1 )B<(t 1 ,t). InEq. @, f L (t-h) and f R (t-h) 
are the convolutions of the density of states function with 
the Fermi-Dirac distributions of the left and right con- 
tacts, respectively, while h(t — ti) is the Fourier trans- 
form of the density of states^ 2 -. The conductance G is 
given by the current divided by the bias voltage V. The 
subsequent time dependent conductance results are com- 
puted by Eq. (U). We will be referring to -q = --^ where 
Ftot = Tl + Tfl, as the asymmetry factor. 

The Kondo effect is a many-body resonance arising 
at low temperatures due to a spin singlet formed from 
the hybridization of the net free spin inside the dot with 
the continuum electrons in the contacts. It originates 
from the seminal work of Jun Kondc 24 , where he discov- 
ered that a divergence in the perturbation series of tun- 
nel couplings yields resistance enhancements in metals 
containing magnetic impurities. Its hallmark is a sharp 
resonance pinned to the Fermi levels of the contacts in 
the dot density of states. The broadening of the Kondo 
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FIG. 2: The grey, red and orange curves (from top to bot- 
tom) show the instantaneous conductance versus time in the 
Kondo timescale after the dot level has been switched to its 
final position for asymmetry factors of 0.85, 0.9 and 0.95, re- 
spectively, at T = 0.009r to t and V = Tk for constant T to t- 

resonance is given by an energy scale Tk (Kondo tem- 
perature), defined as 

Here, D is an energy cutoff, equal to half the bandwidth 
of the conduction electrons, and T to t is the value of the 
total coupling between the dot and the contacts r tD j(e) 
at e = €f- 

In the following, we investigate the instantaneous con- 
ductance for a system in which the dot level is abruptly 
shifted from e\ = — 5r tot to = — 2r tot where T tot —0.8 
eV at t = by a gate voltage, resulting in a transition 
from a non-Kondo state (Tk C T) to a Kondo state. In 
the final state, we infer T K = 0.0025r tot from Eq. ©. 
In the initial short timescale associated with charge fluc- 
tuations, the conductance reaches a maximum for large 
asymmetry factors and then starts to decay, confirming 
previous studies^. Fig. [5] shows the behaviour of the 
instantaneous conductance in the long Kondo timescale 
after the dot level has been switched to its final posi- 
tion. The instantaneous conductance exhibits a complex 
ringing behaviour and the amplitude of the oscillations 
diminishes with decreasing asymmetry factor, completely 
disappearing for symmetric coupling. This is simply be- 
cause the interference between the left contact and the 
Kondo resonance is out of phase with the interference 
between the right contact and the Kondo resonance due 
to opposite signs in Eq. ((4]). The amplitudes of these two 
interference processes are equal for symmetric coupling 
and the oscillations cancel out. 

The frequencies taking part in the oscillations can be 
extracted by a Fourier transform of the time dependent 
conductance. It turns out that there exist two distinct 
frequencies, oj% and u>2, which give rise to a beating be- 
haviour with envelope and carrier frequencies of u\ — 0J2 
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FIG. 3: The red, blue and black curves (from top to bottom) 
represent the instantaneous conductance versus time in the 
Kondo timescale after the dot level has been switched to its 
final position for an asymmetry factor of 0.9 at T — 0.009T to t 
when the bias is equal to V = Tk, V = 4Tk and V = STk, 
respectively, for constant Ttot- 



and ui + ll>2, respectively. We find 0J\ = 2.15cj2, which 
corresponds to the ratio of the distances of the peaks at 
-2.08 eV and -4.46 eV in Fig. [T] from the Fermi level. 
This result demonstrates that van-Hove singularities can 
be probed by the transient current in the Kondo regime. 

The effect of finite bias on the instantaneous conduc- 
tance is depicted in Fig. [3l The influence of finite bias 
is twofold. First, it reduces the amplitude of the oscilla- 
tions. Second, the decay rate of the oscillations, defined 
as the inverse of the time it takes the oscillations to die 
out, increases, indicating that the oscillations are inher- 
ently related to the formation of the Kondo resonance. 
It should be noted that the split Kondo peak oscillations 
arising at finite bias have no discernible effect on these 
results since their frequency is two orders of magnitude 
smaller. 

Finally, the effect of ambient temperature is displayed 
in Fig. |4] Lowering the ambient temperature influences 
the time dependent conductance results in the following 
ways. First, the decay rate of the oscillations decreases. 
Second, the amplitude of the oscillations starts increas- 
ing, but saturates when the temperature approaches the 
Kondo temperature Tk- Lowering the temperature be- 
low Tk does not alter the amplitude anymore. 

Based on these observations and the fact that disconti- 
nuities in the density of states of the contacts can induce 
interference with the Kondo resonance for an asymmet- 
rically coupled system^, we suggest that the beating be- 
haviour in the time dependent conductance is a result 
of interference between the emerging Kondo resonance 
at the Fermi level and the sharp features in the den- 
sity of states of gold located at —2.08 eV and —4.46 eV. 
Since the electrons in the contacts are assumed to be non- 
interacting, the density of states of the contacts is time 
independent. It always looks like Fig. [TJ As a result, 
sharp features in it are also static. On the other hand, 
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the Kondo resonance is part of the density of states of 
the dot. Consequently, its formation is time-dependent 
and dynamical. This naturally leads to persistence of 
the beating behaviour until the Kondo resonance is fully 
formed. The ratio of u>\ to u>2 supports this scenario. 
Moreover, the reduction of the oscillation amplitudes 
with increasing source-drain bias is due to the fact that 
the Kondo resonance starts getting destroyed and the in- 
terference vanishes. The saturation of the oscillation am- 
plitudes below Tff is a further testimony to the proposed 
mechanism, as the Kondo resonance is best developed 
below this scale. Hence, the interference strength stays 
the same. We want to note that the validity of our re- 
sults in transient regime can be checked by considering a 
second system which has a smaller Tk- This can be ac- 
complished by taking a slightly lower final dot level Cdot- 
The time dependent conductance curves overlap for both 
systems as a function of Tk t. This proves that our calcu- 
lations obey universality which is a hallmark of Anderson 
model. 

In conclusion, we have studied the transient current 
in a single electron transistor consisting of gold con- 
tacts in response to an abrupt switching of its dot level. 
We have used the density of states of gold, as obtained 
from ab-initio electronic structure calculations, as input 
in a many-body calculation. We find that an asymmet- 
rically coupled system exhibits complex oscillations in 
the Kondo timescale. The two distinct frequencies that 
give rise to the observed beating behaviour are found 
to be proportional to the separation between the Fermi 
level and the two sharp features in the density of states. 
We interpret this as an interference between the emerg- 
ing Kondo resonance and the van-Hove singularities. We 
note that more frequencies may mix with the two dom- 
inant frequencies in an actual experiment due to addi- 
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FIG. 4: The green, black, blue, red and purple curves (from 
top to bottom) represent the instantaneous conductance ver- 
sus time in the Kondo timescale after the dot level has been 
switched to its final position for an asymmetry factor of 0.9 
at T = 0.0015r tot , T = 0.0030r tot , T = 0.0060r tot , T = 
0.0090r tot and T — 0.0150T to t, respectively, for constant Y to t 
with a bias of V = Tk- 

tional features in the gold density of states, see Fig. [T] 
However, their effects on the reported findings would not 
be discernible as the amplitudes are negligible as com- 
pared to those captured by our fitting. 



Finally, we feel that the situation discussed in this let- 
ter can be realized with today's experimental technology, 
since ultrafast pump-probe experiments can capture the 
picosecond timescale^. This is close to the time needed 
for forming the Kondo resonance^. It would be greatly 
desirable to confirm our predictions experimentally. 
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